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We present a comprehensive analysis for the determination of the confusion levels 
for the current and the next generation of far-infrared surveys assuming three dif- 
ferent cosmological evolutionary scenarios. We include an extensive model for diffuse 
emission from infrared cirrus in order to derive absolute sensitivity levels taking into 
account the source confusion noise due to point sources, the sky confusion noise due 
to the diffuse emission, and instrumental noise. We use our derived sensitivities to 
Q . suggest best survey strategies for the current and the future far-infrared space mis- 

^ I sions Spitzer, AKARI (ASTRO-F), Herschel, and SPICA. We discuss whether the 

, theoretical estimates are realistic and the competing necessities of reliability and com- 

■ pleteness. We find the best estimator for the representation of the source confusion and 

produce predictions for the source confusion using far-infrared source count models. 
From these confusion limits considering both source and sky confusions, we obtain the 
' optimal, confusion limited redshift distribution for each mission. Finally, we predict 

H \ the Cosmic Far-Infrared Background (CFIRB) which includes information about the 

number and distribution of the contributing sources. 

Key words: cosmology: observations - infrared: galaxies - galaxies: evolution - ISM: 
structure - methods: observational 



1 INTRODUCTION 

The detection limits of space infrared (IR) telescope sys- 
tems are governed by instrumental, photon, and confusion 
noise sources. The instrumental noise is intrinsic to the de- 
tector and electronics system, e.g., readout noise and dark 
current fluctuations. The photon noise is due to Poissonian 
fluctuations in the photon counts from the sky background 
emission and the telescope thermal emission which can be re- 
duced by decreasing the telescope temperature. On the other 
hand, the confusion noise is due to both the superposition of 
point sources in crowded fields and the brightness fluctua- 
tion of extended structures at scales of order of the telescope 
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beam size. In the far-IR bands, due to the relatively large 
beam sizes, astronomical observations are mostly affected 
by the sky confusion due to the Galactic cirrus structure 
IIGautier et al. 19921 and the source confusion due to the 
unresolved extragalactic sources. 

Though the confusion level is the unavoidable limit 
for the detection of faint point sources, the source struc- 
ture below the confusion limit also creates background 
fluctuation in the observed image. This means that the 
source distribution can be obtained by measuring the 
background fluctuation which results from a convolu- 
tion of the background sources and the telescope optics 
ULagache fc"Tuget 2000b Matsu hara et ai. 20 00). After the 
successful Infrared Astronomy Satellite mission (IRAS) 
jSoiter, Houck fc JNeugebauer 1987| l, subsequent space IR 
telescopes such as the Infrared Space Observatory (ISO) 
IIKessler et al. 199611 and Spitzer ( |GaUagher et al. 2002| 
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IWerner et ai. 200411 . have increased the levels of sensi- 
tivity and spatial resolution at far-infrared wavelengths. 
In addition, the Japanese AKARI (formerly known as 
ASTRO-F) satellite was launched successfully on Februrary 
21st 2006 jM urakami 1998 Nakagawa 2001 Shibai 2004 
IPearson et al. 2004. Pearson et al. 2006bil , and will perform 
an all sky survey in 4 far-IR bands with much improved 
sensitivities, spatial resolutions and wider wavelength cov- 
erage than the previous IRAS survey over two decades 
ago. Moreover, within the next decade, the Herschel Space 
Observatory (Herschel) HPilbratt 20021 IHarwit 2004ll and 
the Space Infrared Telescope for Cosmology and Astro- 
physics (SPICA) fNakagawa et al. 2004| l will observe the 
much deeper universe with large aperture sizes of 3.5 m. In 
our previous paper (Jeong et al. 2005, hereafter Paper I), we 
numerically generated a high spatial resolution map of the 
Galactic cirrus and investigated the effect of sky confusion 
for current and future space far-IR missions. In addition to 
the simulated sky confusion noise, this work includes real- 
istic source confusion effects assuming various cosmological 
evolutionary scenarios. We also obtain the expected Cosmic 
Far-IR Background (CFIRB). 

The fluctuations in the surface brightness of extended 
structure on scales smaller than the resolution of the 
telescope/instrument beam can produce spurious events 
that can be easily mistaken for genuine point sources, 
since the existence of a source is usually simply derived 
from the difference in signal between the on source and 
some background (off source) position. Such extended 
structure is observed in wispy neutral interstellar dust in 
the Milky Way that is heated by the interstellar radiation 
field and is known as the infrared cirrus lILow et al. 1984|l . 
and is one of the main noise sources in the far-IR 
range. Cirrus emission peaks at far-IR wavelengths but 
was detected in all 4 IRAS bands. There have been 
realistic estimations of the sky confusion from obser- 
vational data by IRAS and ISO fGautier et al. 19921 
IHelou fc Beichman 19*901 [fferbstmcicr ot al. 1998 

[Kiss ot al. "20011 [Kiss, Klaas fc Lemke 2005| 

[Jcon g ct al. 2005J . The intensity of Galactic cirrus is 
a function of Galactic latitude and is serious for wave- 
lengths longer than 60 /xm. Using high resolution maps 
extended from currently available Galactic emission maps, 
we have estimated the sky confusion noise for various 
missions, based on fluctuation analysis and detailed pho- 
tometry over realistically simulated images. In Paper I, 
we concluded that the sky confusion is expected to be 
two orders of magnitude lower for the next generation of 
space missions with larger aperture sizes such as Herschel 
and SPICA, but on the other hand, current 60-90 cm 
aperture missions such as Spitzer and AKARI, will have 
to endure a sky dominated by confusion noise at long 
wavelengths (> 100 ^m). In this paper, we also discuss 
the composite effect of sky confusion and source confusion 
from the implementation of a realistic model for the source 
distribution. 

Note that additional components affecting the confu- 
sion noise are the presence of faint, unresolved, asteroids 
and the zodiacal light. Though bright asteroids can be easily 
rejected via source confirmation between 2 pointings/scans 
over the same sky position, faint asteroids may contribute 
to fluctuations assumed to be from faint galaxies. Although 



Abraham et al. (I1997II found no such fluctuations on small 
scales from the ISOPHOT (photometer on board the ISO) 
data, more sensitive instruments may detect them. Thus, 
these two components may show a non-negligible contribu- 
tion to the far-IR confusion near the ecliptic plane. 

This paper is structured as follows. In Section ^ we 
briefly describe the confusion noise due to sky brightness 
fluctuations and the source confusion due to extragalactic 
point sources. We present our source count models includ- 
ing galaxy evolution and the simulated images in Section 
13 Based upon the specifications of each IR mission, we es- 
timate the source confusion noise through fluctuation and 
photometric analyses in Section 2] We present the expected 
results from our simulations in Section 13 Our conclusions 
are summarised in Section |S| 

Throughout this work we assume a concordance cos- 
mology of Ho = 72 km s"^Mpc"\ = 0.3, and Q.a =0.7, 
unless otherwise explicitly stated. 



2 CONFUSION DUE TO EXTRAGALACTIC 
SOURCES 

The galaxy confusion limit is defined as the threshold of the 
fluctuations in the background sky brightness below which 
sources cannot be discretely detected in the telescope beam 
9 ^ A/D, where A is the observation wavelength and D is the 
telescope diameter. Thus, the fluctuation noise arises from 
the same origin as the galaxies that one is aiming to detect. 
If we assume galaxies are distributed as a power law in flux, 
5", down to some detection limit Sum, 

iV(>5) = iVH„,(^)"", (1) 

where N{> S) is the number of sources per unit solid angle 
with flux greater than S, a is the slope of the integral source 
counts (where a = 1.5 for a Euclidean Universe) and Num 
is the number of sources brighter than Sum per unit solid 
angle. Assuming that the counts flatten at some faint flux, 
5*0, i.e. a{So) ~ 0, then the intensity of the background (in 
Jy/sr) up to some maximum flux, S'max, corresponding to 
these sources is given by. 

The fluctuations contributed by sources below the detection 
limit Slim is given by the second moment of the differential 
source counts dN/dS, osc. 

Assuming the power law distribution of sources given in 
equation Q, equation ^ can be evaluated to give 




For the Euclidean case, the dominant sources contributing to 
the background intensity are those just below the detection 
limit Slim IIMatsuhara et al. 2000l|Lagache fc Puget 2006bl l. 
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However, the strong evolution detected in the galaxy pop- 
ulation steepens the source counts and produces super Eu- 
clidean slopes (q > 1.5) and the sources around the detec- 
tion limit also contribute significantly to the fluctuations in 
the background intensity. 

Rigorous theoretical definitions of confusion have been 
presented by Scheuer (I1957II and Condon (1974J. Hogg 
1120011 has highlighted more practical aspects of galaxy 
confusion noise. An analytical derivation broadly fol- 
lowing Franceschini et al. 1)1989^ is given below. Note 
that the clustering of sources will complicate the con- 
fusion noise [e.g., as in the case of radio sources, 
Condon H1974|l ] although here, for clarity, we do not 
treat this effect. Some authors lIFranceschini et al. 19891 
ITakeuchi fc Ishii 20041 [Negrello et al. 2004| | have investi- 
gated the effect of the clustering of sources on the confu- 
sion limit. They found that clustering possibly increases the 
level of the source confusion by 10% for AKARI mission 
l|Takeuchi fc Ishii 2004^ and bv 10 ~ 1 5% for Spitzer and 
Herschel missions l |Negrello et al. 2004| l in the far-IR range. 

Assuming that the sources are distributed randomly 
over the sky described by a power law form N{S) oc and 
a corresponding differential distribution given by dN/dS = 
kS ' where 7 = a — 1, the detector response to a source 
of flux S, at a position [6, (j)) from the axis of a detector 
beam of proflle (point spread function) f{6, (j)) is given by 
X = Sf{6,(f)). Hence the mean number of responses with 
amplitudes between x, x + dx in a. solid angle dQ, is given 
by; 



R{x) 



n{x/f{e, 4>)) 



dQ, 



(5) 



where fli, is the solid angle of the beam in steradians. Note 
that for the power law distribution of sources discussed 
above, equation ^ can be rewritten as, 



R{x) = kx 



f{9,(f,y'^dn = kx-^'Qe 



(6) 



where f2e = f f{9)~'~^dQ, is the effective beam 
size llCondon 19741 . Taking the second moment of the R{x) 
distribution (the variance) gives the fluctuation of the re- 
sponse, (Tsc ; 



X R{x)dx, 



(7) 



where Xc is a cut off response introduced to stop the variance 
from diverging at bright source fluxes. More practically the 
confusion limit Xc (corresponding to a cut off flux Sc) is set 
to some factor of the confusion noise such that Xc ~ qo'sc, 
where the factor q (usually chosen values between 3 and 
5) limits the largest response included in the evaluation of 
the confusion noise CTsc [values of ~ 5 are assumed in the 
calculations of Franceschini et al. 1199111 ] . The difference in 
assuming a cut off in the response as opposed to a cut off 
in flux is that weak contributions from strong sources are 
not neglected, as even a strong source far from the axis of 
the beam may contribute significantly to the point spread 
function of the beam. 

Assuming for clarity in the calculations, a circular Gaus- 



-4(9/60) ln2 



such 



sian beam profile, /(e) = f{{9/eo) ) 
that dQ = 29 d9/9l where 9o is the Full Width at Half Max- 



imum (FWHM) of the beam, integrating equation (|7|l over 
the solid angle of the beam gives; 



cr^{xc) = Tv9aP{xc) 



nP{Xc 



(8) 



(9) 



where P, effectively the power in the fluctuations, is given 
by, 



P{Xc 



x'dx n(x//(e))e''®'"''de 



(10) 



Thus, the confusion limit can be directly related to the 
FWHM of the instrumental beam. For the simple power law 
representation of the distribution of extragalactic sources 
given previously and the definitions of equations Q and 
JBJ, the confusion limit is given by. 



J3-7)/2 



- 7 



(11) 



Therefore, the confusion noise limit will be a complex func- 
tion of the beam size 9o, the source counts N{S), the cut off 
in flux Sc or response Xc and the factor q. For the assumed 
symmetric Gaussian beam profile, CTc tx 9o^'^~' 

The confusion due to faint galaxies will be worse at 
longer wavelengths and smaller telescope diameters. Since 
the confusion noise is related to the mean number of re- 
sponses (the source density) and the cut off response qjx^, 
a useful, practical benchmark for the confusion limit can 
be set by limiting the number sources per beam before the 
beam becomes confused. Ideally, the confusion limit would 
be determined by the fiux at which the source density be- 
comes 1 source per beam although more realistically a limit 
of between 1/20-1/50 sources per beam (20-50 beams per 
source) is assumed [e.g., Hogg (2001)]. Note that this differ- 
ence is due to the contribution from faint sources undetected 
below a certain limit. 



THE INPUT CATALOGUES AND 
SIMULATED IMAGES 



Many authors JCondon 19741 IFranceschini et al. 1989| 
IHelmTXTB eichman 19*901 IGautier et al. 19921 

IHerbstmeier et al. 19981 IKiss et al. 200T1 IKiss et al. 20031 
[Jeong et al. 2005^ have described the source confusion due 
to extragalactic point sources and the sky confusion due to 
Galactic cirrus as separate issues since most extragalactic 
surveys are limited to low background regions. However, in 
order to cover much larger survey areas, we also naturally 
have to include the medium-to-high background regions. In 
addition, since the cirrus structure covers both large and 
small scales, both of these contributions should be consid- 
ered together. By using various source distributions and a 
realistically simulated cirrus background, we address the 
effects of sky confusion and source confusion simultaneously 
with a more realistic framework. 
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Figure 1. Model input spectral energy distributions used for in- 
put catalogues to tlie simulation. A four component model com- 
prising of a normal galaxy, starburst galaxy, ultraluminous galaxy 
and AGN dust torus are included. The source spectral energy 
distributions are based on the models of Efstathiou et al. (2001), 
Efstathiou & Rowan-Robinson (2002), Rowan-Robinson (1995). 



3.1 Input Catalogues 

The input catalogues are prepared using the models of 
Pearson (200 1| (hereafter GPP). GPP is a far-IR model 
based on the IRAS colours and luminosity function of 
galaxies [see also Pearson & Rowan- Robinson (I1996|l ]. The 
model incorporates a 4 component parameterization of 
galaxy Spectral Energy Distributions (SED) segregated 
by IRAS colours ( Rowan- Robinson & Crawford 1989). 
A normal galaxy population modelled on the cool 
100 /im/60 /xm colours identified with infrared cir- 
rus ijLow et al. 19841 lEfstathiou fc R.owan-R.obinson 2003fl . 
A starburst population based upon the warm 100 
/im/60 pm colours of IRAS galaxies with the 
archetypical starburst galaxy M82 as a template 
SED ( [Efstathiou, Rowan-Robinson fc Siebenmorgen 2000at . 
An ultraluminous galaxy popula- 

tion ijSanders fc Mirabel 1996^ representing the high 
luminosity tail of the IRAS starburst galaxy popu- 
lation and representative of the archetypical ULIG 
ARP220 ( [Efstathiou, Rowan- Robinson fc Siebenmorgen 2000at . 
An AGN (Seyfert 1 & Seyfert 2) population modelled on 
a 3-30 nm dust torus component lIRowan-Robinson 199511 
defined by hot 25 fim colours. Though recent 

Spitzer observations showed silicate dust emissions in 
luminous and low-luminosity AGN IIHao et al. 20051 
[Siebe n morgen et al. 2005, ..Sturm et al. 2005^ . we do not 
include it. The input SEDs are shown in Figure These 
SEDs have been shown to provide good fits to the SEDs 
and colours of sources in the /SO-ELAIS survey at 6.7, 
15, 90 and 170 /im (|R owan-Robinson et al. 20041 and also 
recently for the first results from the Spitzer SWIRE Legacy 
survey ijRowan-R.obinson et al. 2005fl . 

To produce the input source distributions, we calculate 
the total number of sources per unit solid angle at observa- 
tion wavelength, Ao, down to some flux limit S\o', 



where f{z) and e{z) are evolutionary functions in the lu- 
minosity and density of the source population respectively. 
The integration is made over the luminosity function (num- 
ber density of objects as a function of luminosity), <j>{L) and 
the cosmological volume V, enclosed inside a limiting red- 
shift z{L, S) defined as the redshift at which a source of 
luminosity, L, falls below the sensitivity, S{Xo)- 

Luminosity functions, 0(i), are determined from the 
All Sky IRAS PSCz catalogue at 60 /im by Saunders 
et al. 1I2OOOI1 for starburst and normal galaxies follow- 
ing colour criteria akin to those of Rowan-Robinson & 
Crawford C 1 9_89 5 • Similarly, the hot AGN population is 
well represented by the 12 /im sample of Rush et al. 
119931 using the luminosity function of Lawrence et al. 
f 1986"). In addition to the above, following CPP, we 
utilize a log exponential luminosity function, defined at 
60 fim, to represent the ULIG population, [referred to 
as the Burst model] which was originally implemented 
to address the paradigms of the strong evolution in the 
galaxy source counts observed with the ISO in the mid 
to far-IR, at sub-mm wavelengths with the Submillimetre 
Common User Bolometer Array (SCUBA) and the de- 
tection of the CIRB at ~ 170Atm ([Oliver et al. 1997| 

Small, Ivison fc Blain 1997[ iKawara et al. 1998] 

Hughes et al. 1998[ IFlores et al. 19991 Et icri et al. 199i| 
Aussel et al. 19991 [Gruppioni et al. 1999| 

Puget et al. 1999[ lEfstathiou et al. 2000bl 

Biviano et al. 20001 lElbaz et al. 20001 ISerjeant et al. 20001 
Lag ache et al. 2 00qa| [Matsuhara ct al. 20001 

Scott et al. "200211 . This model was found to provide a 
good fit to both the number counts and redshift distribu- 
tions of galaxies from sub-mm to near infrared wavelengths 
as well as the cosmic infrared background. To shift the 
luminosity function from the wavelength at which the 
luminosity function is defined, Alf, to the observation 
wavelength, Ao, the ratio I/(Ao)/-£/(Alf) is obtained via the 
model template spectra. 

We prepare 3 types of input catalogue and produce 
many versions of each catalogue to reduce statistical errors. 

(i) No-Evolution Model - No evolution is assumed for all 
galaxy components. 

(ii) Luminosity Evolution Model - Luminosity evolution is 
included following Pearson & Rowan- Robinson (1996). Only 
evolution in the luminosity of the source population is as- 
sumed in this scenario. This luminosity evolution follows a 
parametric form as a function of redshift of f{z) — {1 + z)^ , 
where k is independently defined for each galaxy type to 
produce the best fits to the far-IR counts. 

(ill) Burst Evolution Model - Luminosity and density evo- 
lution is included following the CPP model. The assumed 
evolution in both luminosity and density incorporates a 
burst in a specific redshift range and follows a paramet- 
ric form as a function of redshift of f{z) and e(z) = 1-1- 
kexp [— (z — Zp)^/2(.<j^] , where the parameters k, Zp, ui are 
independently defined for the luminosity and density evolu- 
tion and galaxy type to produce the best fits to the far-IR 
counts. 

Details of the evolution included in the models can be 
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found in CPP. Figure |5| shows the evolutionary source count 
models for the Spitzer far-IR galaxy counts at 70 nm and 
160 fj,m respectively ijDole et al. 2004a^ . The two evolution- 
ary models provide reasonable fits and also cover the spread 
in uncertainty shown in the observations. Therefore, in this 
work we consider the burst evolution and luminosity evo- 
lution models as upper and lower limits to the numbers of 
far-IR sources respectively. We do not attempt to fit the 
mid-near-infrared Spitzer counts. For specific updated evo- 
lutionary models, please see Pearson Il2005t and Pearson 
Il2006all . 

3.2 Simulated Images 

Based upon the 3 types of input catalogue, we generated 
simulated images in each band for various space missions. 
The image size for the distributed source simulation is ~ 

5.3 square degrees. In order to check the effect of sky con- 
fusion noise, we also include high resolution cirrus maps 
whose mean cirrus brightness ranges from 0.3 MJy sr~^ to 
25 MJy sr~^ at 100 ^m as the background by using the 
method described in Paper I. 



4 SIMULATION RESULTS 

As we described in Section |5| there are many definitions 
for the source confusion. We compare these definitions and 
propose an appropriate definition for source confusion. In 
Table0 we list the basic instrumental parameters of present 
and future IR space missions; the aperture of the telescope, 
FuU Width at Half Maximum (FWHM) of the beam profile 
and the pixel size for each detector. 

4.1 Definition by 'Beams per Source' 

First, we estimate the source confusion from the classical 
definition of source confusion, 'beams per source' (sources 
per beam). Though often cited in the literature as 'sources 
per beam', we use the term 'beams per source' for our defi- 
nition of source confusion in order to compare the effect of 
source confusion according to the different beam sizes for 
each mission. We check the source confusion by changing 
the number of beams per source criterion from 10 to 50. 
Fig. 13 shows the source confusion assuming this definition 
for each mission and each model. Hogg f2001^ showed that 
30 beams per source is the minimum photometric criterion 
where the source counts are steep, and suggested 50 beams 
per source for the definition of source confusion. Rowan- 
Robinson 1200111 adopted 40 beams per source. In recent 
papers. Dole et al. (2003, 2004b) suggested that source con- 
fusion could be defined by a source density criteria, cor- 
responding to ~ 12 beams per source for the Spitzer mis- 
sion. They have estimated a source confusion limit for the 
Spitzer mission of 3.2 mjy for SW band and 40 mjy for 
LW band, respectively. These results are similar to our esti- 
mations with the definition of 12 beams per source for the 
luminosity evolution model in the SW band and for the burst 
evolution model in the LW band, respectively. These consis- 
tences suggest that the source distribution model used for 
the estimation of the source confusion limit for the case of 



Spitzer ( |Lagache, Dole fc Puget 2003| [Lagache et al. 2004t 
also predicts that a starburst component, and ULIG com- 
ponent dominate the SW and LW bands respectively. In 
Table|5| we list the source confusion limits estimated by the 
definition of 12 and 40 beams per source. As seen from Fig. 
13 we found that the source confusion did not increase at 
a constant rate and that the slope of the source confusion 
slightly varies according to the source distribution model 
and the resolution of the mission, especially for burst evo- 
lution model. Therefore, we conclude that we can not apply 
the same definition of 'beams per source' in a generic way 
for different missions. 

4.2 Definition by Fluctuation 

Another criterion for the quantification of source 
confusion can be defined by the fiuctuation from 
beam to beam due to the distribution of the 
point sources (lUondon 19741 ll^'ranceschini et al. 1989| 
Vai sanen et al. 20011 IXu et al. 200111 . Since the beam size 
is large and the source counts are steep, the usual definition 
of 'beams per source' may not be valid in the case of far-IR 
photometry. For the PSF, we assume an ideal circular 
aperture Airy pattern corresponding to the aperture size 
of the telescopes for each mission except for AKARI for 
which we use the theoretical PSFs estimated from the tele- 
scope design Peong et al. 200'3| . To visualize the iteration 
procedure, we plot the 5'iim/o'sc (5'iim) ratio as a function of 
Slim for the Spitzer and Herschel & SPICA missions for the 
case of the burst evolution model (see Figure For the 
SW band of the Herschel & SPCIA missions, the Sum/fsc 
ratio is always greater than 5, which means that we can 
not obtain a solution for source confusion, even for g = 5. 
Dole et al. 120031 estimated the source confusion with their 
best estimator for the Spitzer mission. With optimized q 
parameters of 3.8 in SW band and 6.8 in LW band, they 
obtained source confusion limits of 3.2 mJy and 36 mJy at 
70 /im and 160 /xm, respectively. Using our source count 
models and the same q parameters, we obtained the source 
confusion limits of 3.7 mJy and 21 mJy for the luminosity 
evolution model and 12 mJy and 60 mJy for the burst 
evolution model at the same wavelengths. As we showed in 
Section l4.ll though the source density of the source count 
model used in Dole et al. 12003t is similar to the luminosity 
evolution model in the SW band and the burst evolution 
model in the LW band, we find that the fiuctuation due to 
distributed sources in the burst evolution model is much 
stronger than that of the model of Dole et al. f2T)033. This 
amount of fluctuation may degrade the flnal detection 
limits. 

In Table 13 we list the source confusion limits estimated 
from the definition by fluctuation for the cases of 5 = 3 and 
g = 5. For q = 5, there are no solutions for the SW band of 
Herschel & SPICA missions. Even in the LW band, we could 
not find a reasonable solution since the estimated source con- 
fusion for the two evolutionary models gave identical results. 
However, when we attempted photometry on the simulated 
images including the point sources for each source count 
model, we encountered some limits to the source detection. 
Therefore, we conclude that we can not apply a constant 
q to the estimation of source confusion generically for all 
cases. 
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Table 1. Instrumental parameters for various space missions. 



Aperture Central Wavelength FWHM " Pixel size 
(meter) (arcsec) (arcsec) 



Jpace Mission 




SW 


LW 


SW 


LW 


SW 


LW 


Spitzer ^ 


0.85 


70 


160 


17 


35 


10 


16 


AKARI = 


0.69 


75 


140 


23 


44 


27 


44 


Herschel 


3.5 


70 


175 = 


4.3 


9.7 


3.2 


6.4 


SPICA 


•■^.5 


70 


160 


1..'^ 


9.7 


1.8 


.■^.6 



FWHM of diffraction pattern. 

^ Two MIPS bands for the Spitzer mission have 3 bands with central wavelengths of 24 //m, 70 /im and 160 /im. We use 70 //m band 
SW band and 160 //m band as LW band in this paper. 

= AKARI/FIS (Far Infrared Surveyor) has a WIDE-S band in SW and a WIDE-L band in LW. 

PACS has a 'blue' array in two short wavelength bands (centered at 70 fim and 110 /im) and a 'red' array at longer wavelengths 
(centered at 175/im). In this paper, we use only the 70 ixm. band of the 'blue' array. 

'' '' Since one of our motivations in this paper is to compare the confusion limits, we use only common bands among all far-IR bands 

the considered space missions. 

Note that though the central wavelength of the 'red' array in PACS is 175 ^m, our estimated results arc for 160 /im. Due to this 
wider beam, the confusion noise is expected to increase by ~ 15% at 175 jxm compared to that at 160 jum. 



Table 2. Source confusion estimated by the definition of 12 and 40 beams per source. 



No evolution 
(mJy) 

SW LW 
Space Mission 12b 40b 12b 40b 



Spitzer 0.53 
AKARI ^ 1.9 
Herschel & SPICA 0.014 



1.9 


8.8 


23 


6.0 


13 


34 


0.066 


0.65 


2.1 



Luminosity evolution 
(mJy) 

SW LW 



12b 


40b 


12b 


40b 


2.1 


5.1 


17 


34 


5.1 


12 


23 


47 


0.10 


0.47 


2.6 


6.8 



Burst evolution 

(mJy) 
SW LW 



12b 


40b 


12b 


40b 


6.9 


13 


42 


69 


14 


26 


52 


87 


0.11 


0.65 


4.3 


20 



" source confusion defined by Dole et al. (2004b): flux corresponding to 12 beams per source. 

source confusion defined by Rowan-Robinson (2001): flux corresponding to 40 beams per source. 
= Wide-S band for SW and Wide-L band for LW. 



Table 3. Source confusion estimated by the definition of fluctuation with q = 3 and q = 5. 



No evolution Luminosity evolution Burst evolution 

(mJy) (mJy) (mJy) 

SW LW SW LW SW LW 



Space Mission 


q = 3 


q = 5 


q = 3 


q = 5 


q=3 5=5 


q = 3 


q = 5 


9 = 3 


q = 5 


q = 3 


q = i 


Spitzer 


0.065 


0.28 


3.8 


10 


0.51 1.9 


14 


33 


0.51 


1.9 


41 


87 


AKARI 


0.36 


1.4 


5.6 


15 


2.4 6.8 


20 


43 


2.4 


23 


59 


115 


Herschel & SPICA 


a 


0.002 


0.042 


0.34 


0.006 


0.18 


1.5 






0.18 


1.5 



" no solution for this case. 
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Figure 2. Comparison between source count models and observations for SW (left) and LW (right) band. All observations except for 
FIRBACK field (Dole et al. 2001) were obtained from Multiband Imaging Photometer (MIPS) for Spitzer (Rieke et al. 2004) in the 
Chandra Deep Field South (CDF-S), the Bootes field and Marano field (Dole et al. 2004a). The observational mode (scan map) provides 
multiple sightings of each source, typically 10 and 60 at 70 ^m in the Bootes and CDF-S, respectively. Because of the low redundancy level 
of the 160 ^m data in the Bootes field, they exclude those data. Upper and lower lines show our burst evolution model and luminosity 
evolution model, respectively, with the hashed area defining upper and lower bounds of far-IR source distributions assumed in this work. 
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Figure 3. Source confusion with the definition of beams per sources for AKARI, Spitzer, and Herschel & SPICA missions. Each line 
shows the source distribution model used in this paper. The two vertical lines show 12 beams (left) and 40 beams (right) per source, 
respectively. 



4.3 Composite Definition by Fluctuation and 
Photometry 

In an attempt to define the source confusion for all cases, 
we implement a composite definition by fluctuation and pho- 
tometry. As a first step, we attempt photometry on the sim- 
ulated images created from each source count model. The 
definition of completeness and reliability for quantifying the 
source detection efficiency is widely used in photometry. The 
'completeness' is defined as the fraction of detected sources 
to the original input catalogue sources and the 'reliability' 
as the fraction of real sources to all detected sources. For our 
photometry, we removed the false or spurious detections by 



comparing the positions of detected sources with those in 
the input catalogue. The real sources as defined in the reli- 
ability, refers to sources of which the measured flux agrees 
with the input flux to within a 20% error. Note that the 
reliability at a given flux range is not always higher than 
the completeness. An excess of sources near the detection 
limit, or more likely an overestimation of the flux of sources 
at or near the detection limit can be caused by a step effect 
where the underlying, unresolved sources are entering the 
PSF, affecting the sky-subtraction. We carried out aperture 
photometry on the simulated images using the SExtractor 
software v2.2.2 IIBertin fc Arnouts 199611 . The most influ- 
ential parameters are the size of background mesh and the 
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Table 4. Source confusion estimated by the composite definition 
of both fluctuation and photometry. 
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Figure 4. Sw^/osc ratio as a function of for Spitzer (upper) 
and Herschel &i SPICA (lower) missions in the case of a burst 
evolution model. We also plot Sh^/cTsc = 3 (dashed line) and 
5 (dashed dotted line). In the case of the Herschel & SPICA 
missions in the SW band, note that Siim/csc is always greater 
than 5, and the plateaus in the SW and LW band arc due to the 
bump seen in the source counts. 



threshold for the source detection for the aperture photome- 
try. Since the size of background mesh is related to the range 
of scales of the background fluctuation, as we set a smaller 
mesh size, we can detect smaller fluctuation. Therefore, if 
we use smaller size of background mesh and thresholds, we 
obtained many more spurious sources. We set the size of 
the background mesh to ~ 2 times the measuring aperture, 
and the detection threshold as 4, which is optimized for bet- 
ter reliability of the detected sources thus reducing the false 
detection rate. 

As discussed in Section r4.2l we could not use a constant 
value of q for the estimation of source confusion. Current 
and future space missions, will detect much fainter sources 
to higher sensitivity, which means that we will observe high 
source densities to extremely faint detectable flux levels. 
Therefore, a significant factor contributing to the source de- 
tection comes from both the faint sources below the detec- 
tion limit and the high source density above or around the 
detection limit. With this assumption, we include the con- 
tributed fluctuations from the sources above the detection 
limit as well. In order to find the limiting fiux affecting the 
source detection, we use a photometric method on the sim- 
ulated images. We set the limiting fiux Sum in equation |21 
to be the flux that the completeness reaches ~ 80% where 
the completeness level abruptly falls off. We do not use the 
reliability criterion in this deflnition since the reliability can 
be improved by using optimized photometric methods. We 
also assume that the sources above this flux level do not con- 
tribute signiflcantly to the source confusion. We obtain the 
flnal source confusion from a 4a" fluctuation in order to com- 
pare with the threshold used in the photometry. In Table |11 
we list the source confusion limits from our best estimators. 

Though we could not obtain the results for the SW band 
of Herschel & SPICA missions with the original deflnition of 
the fluctuation with constant q, we could obtain the results 
by our new best estimator. We also find that the source con- 
fusion limits by our definition are mostly consistent with a 
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Burst evol. 




(m 


Jy) 


(m 


Jy) 


(mJy) 


Space Mission 


SW 


LW 


SW 


LW 


SW LW 


Spitzer 


0.48 


9.4 


2.0 


23 


6.5 66 


AKARI 


2.0 


14 


6.3 
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17 87 


Herschel & SPICA 


0.0095 


0.63 


0.077 
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Figure 5. Comparison of sky confusion between Kiss et al. (2005) 
and Jeong et al. (2005) for the AKARI mission at 170 /^m. Two re- 
sults shows significant differences in the range of low and medium 
cirrus brightness. 

completeness of ~ 65% for all missions. Therefore, we con- 
clude that our definition can explain the behaviour of source 
confusion well in the far-IR range, regardless of individual 
space mission characteristics. 



4.4 Predicted Confusion Limits for Current and 
Future Missions 

In our study of sky confusion in Paper I, we predicted the 
confusion limits considering both sky confusion due to cir- 
rus structure and source confusion for each mission from a 
rather simple (power law) source distribution models and 
plotted the sky confusion levels for an assumed range of av- 
erage cirrus brightness {B\) in Table |^ (see details in Paper 
I). Kiss et al. i2005J have also estimated the all-sky cirrus 
confusion for various IR missions. Fig.|5]shows a comparison 
of sky confusion between Kiss et al. llSOtlSTi and Jeong et al. 
(2005 ) for the AKARI mission at 170 /im. Since the Kiss et 
al. (2005) model assumes a linear relationship between the 
cirrus brightness and sky confusion, it often gives different 
results in low to medium cirrus brightness. 

Since, in general, the cirrus fluctuations are not repre- 
sented by Gaussian noise, we can not directly sum the two 
noise contributions (sky confusion and source confusion) for 
the estimation of the flnal confusion noise. However, in a 
recent paper, Heraudeau et al. 112004 ll found that the distri- 
bution of cirrus fluctuation is near-Gaussian, at least at the 
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Table 5. Estimated confusion limits due to Galactic cirrus for 
each mission. 



Sky Confusion Limits (mjy) for (-Ba)" 





0.5 MJy/sr 


5 MJy/sr 


15 MJy/sr 


Space Mission 


SW 


LW 


SW 


LW 


SW LW 


Spitzer 


0.10 


7.1 


1.2 


10 


11 19 


AKARI 


0.73 


9.4 


5.4 


18 


27 42 


Herschel &c SPICA 


0.003 


0.16 


0.027 


0.21 


0.31 0.37 



° mean brightness of cirrus emission 

ISOPHOT 90 and 170 fim spatial frequencies. Thus, we use 
equation 1131 to estimate final confusion limits 5conf in this 
paper by summation of the two confusion noise components: 

5'conf = Sec ~t~ 'S'sc 5 (i-"^) 

where Sec is the sky confusion limit due to cirrus structure 
and Sac is the source confusion limit. Fig. [5] shows the final 
confusion limits due to both sky and source confusion for 
each mission and each evolutionary model. The sky confu- 
sion will be negligible for the Herschel and SPICA missions 
even in regions of high cirrus brightness and the dominant 
confusion contribution is predicted to be the source con- 
fusion in high Galactic latitude regions for the purpose of 
cosmological studies. However, care must be taken when con- 
sidering the Spitzer and AKARI missions in the mean cirrus 
brightness range (-Ba) > 20 MJy sr^^ for the LW band. In or- 
der to check the effect caused by the combination of both sky 
and source confusion, we attempt photometry on the simu- 
lated images including various cirrus background levels and 
source count models and compare them with the photomet- 
ric results considering only source confusion. These results 
are presented in Appendix IXI 

We summarize the confusion limits including both sky 
and source confusion for each mission in Table M for the 
low cirrus background regions dominated by source confu- 
sion. The mean brightness used in the estimation of sky 
confusion is 0.5 MJy sr"^ and 1.0 MJy sr"^ for the SW 
and LW bands, respectively. The sensitivity for the AKARI 
mission is estimated based upon the recent laboratory ex- 
periments, and various detector characteristics and obser- 
vational environments which may affect the final sensitivity 
(Jeong et al. 2003, 2004; Matsuura et al. 2002; Shirahata et 
al. 2004). Though the single scan sensitivity in survey mode 
of the AKARI mission is much lower than the confusion 
limit, the sensitivity at higher ecliptic latitudes where many 
scans overlap should approach the confusion limit. For the 
Herschel mission, since the sensitivity is 3 mJy, we expect 
that most observations in the SW band will be dominated 
by instrumental noise and those in the LW band will be 
near the confusion limit. However, the SPICA mission will 
achieve the confusion limit in both bands since it will have 
a large aperture telescope cooled to very low temperatures 
jNakagawa et al. 2004| i. 

The detection limits as a function of integration time 
can be an important indicator of the confusion level. Based 



upon recent hardware measurements for the AKARI mission 
and the Spitzer results from Dole et al. (2004b), we compare 
the detection limits as a function of an integration time in 
Figure [7| In case of the AKARI mission, we also plot the 
results of different sampling modes for the high background 
regions. The decreasing component of the noise is the usual 
behavior of noise with increasing integration time while the 
flattening component is mainly due to the confusion. Since 
a low cirrus region was assumed in this estimation, the flat- 
tening of the noise results from the source confusion. In the 
SW band, the detection limits are mostly determined by the 
instrumental noise and the source confusion noise irrespec- 
tive of the amount of cirrus background except for the very 
highest cirrus regions, i.e. near the Galactic plane. How- 
ever, in the LW band, we find that the sky confusion affects 
the noise level severely in the regions brighter than ~ 20 
MJy sr"^ for the AKARI mission and ~ 30 MJy sr"^ for 
the Spitzer mission, respectively (see also the plot of con- 
fusion limits in Fig. [SJ. If the different sampling mode for 
high background regions available for the AKARI mission is 
used, the detection limits are mostly limited by instrumen- 
tal and photon noises in the SW band unless the exposure 
time is very long. However, in the LW band, we found that 
the sky confusion limits compete with the instrumental and 
photon noises even for short exposure observations. 



5 EXPECTED RESULTS 

Based upon two evolutionary (and a no evolution) scenarios 
of source distributions, we have estimated the expected con- 
fusion limits for various IR missions in Section f4.4l The two 
evolutionary models, the luminosity evolution model and 
burst evolution model, provide upper and lower limits for 
the final detection limits. Here, we discuss the results ex- 
pected from the two evolutionary models from our estimated 
sensitivities. 

5.1 Expected Optimal Redshift Distribution 

Once the realistic detection limits have been determined, we 
can obtain the expected optimal confusion limited redshift 
distribution from the models. The number-redshift distri- 
bution to our limiting flux Sum, can be extracted directly 
from our simulated catalogue. In Figure |H] we show the ex- 
pected redshift distributions and expected number of sources 
per square degree for each mission at the predicted con- 
fusion limit. We consider this as the optimal redshift dis- 
tribution for each mission. Since a fairly large number of 
sources below z < 1.0 are detected in the LW bands for the 
no evolution model, the redshift distribution for each mis- 
sion has a significant difference depending on whether there 
is evolution or not. Though we can find the peaks in the 
SW band of the burst evolution model and the LW band 
of luminosity evolution model for the Spitzer and AKARI 
missions, we still have difficulty in finding a significant dif- 
ference between the two evolutionary scenarios. However, 
we expect that the large 3.5m aperture missions (e.g. Her- 
schel and SPICA) will have advantages over the 60-90cm 
aperture missions {Spitzer and AKARI) in distinguishing 
between the evolutionary scenarios due to the larger num- 
bers of sources detected in the redshift range of 1-3. Since 



10 W.-S. Jeong et al. 



10 

1 

^ 10 
03° 1 

0.10 

0.01 

0, 



AKARI 1 




1 


I Burst evol.: 

! -No evol. 


Spitzer i 












Herschel & SPICA | 






@sw - 





1.0 

mean brightness [MJy/sr] 



10.0 



1000 
100 
10 

1 100 

)" 

10 
10 



AKARI 



Burst evol.j 
Lum. evol.: 
No evol. 



Spitzer 



Herschel & SPICA 



@ LW 



10 

mean brightness [MJy/sr] 



100 



Figure 6. Confusion limits considering both sky confusion and source confusion. The two vertical lines show the mean cirrus brightness 
1.0 MJy sr^^ for SW band (left) and 15 MJy sr~^ LW band (right), respectively. In dark cirrus regions, source confusion is expected to 
dominate. Note that sky confusion increases as the mean cirrus brightness becomes larger and the source confusion has a constant value 
irrespective of the mean cirrus brightness. 



Table 6. Final confusion limits considering both source confusion and sky confusion. 



Sensitivity No evolution Luminosity evolution Burst evolution 

(mjy) (mjy) (mjy) (mjy) 



Space Mission 


SW 


LW 


SW 


LW 


SW 


LW 


SW 


LW 


Spitzer 


6 


15 


0.49 


12 


2.0 


25 


6.6 


67 


AKARI " 


12 (200) 


12 (400) 


2.2 


18 


6.4 


34 


17 


88 


SPICA ^ 


3 


3 


0.010 


0.66 


0.077 


2.3 


0.10 


4.4 


Herschel 


3 


3 


0.010 


0.73 


0.077 


2.6 


0.10 


5.1 



° 5(T sensitivity for slow scan mode (scan speed: 8 arcsec/sec, reset: 1 sec) and survey mode (parenthesized values, 1 scan). The 
sensitivity of the survey mode at high ecliptic latitude regions can be improved due to the high redundancy. 

^ 5a sensitivity without source confusion and sky confusion. In the case of the Herschel & SPICA mission, we commonly use the 5(t 
sensitivity of Herschel mission ( Poglits ch et al. 2002^ . 
Our actual estimations in the LW band are for those in 160 /im. We considered the effect of wider beam at 175 ^im of the LW band. 



the source confusion should be more severe for stronger evo- 
lutionary models, the total number of detected sources in 
the no evolution model is largest in the LW band except for 
the Herschel and SPICA missions. 

From Figure |H1 we show that the confusion limit for 
the Spitzer and AKARI missions, restricts the majority of 
the sources to redshifts of around or less than unity. For 
these missions, we could expect of the order of a few tens of 
sources per redshift bin per square degree out to redshifts 
of 2-3 (Note that the MIPS 70 /xm band is more sensitive in 
theory but suffers from detector array problems in practice) . 
However, for the Herschel and SPICA missions, we could ex- 
pect 500- 1000s of sources per square degree per redshift bin 
out to redshifts of 2-3. This would allow direct probing of the 
dusty starburst populations detected at 24 fim with Spitzer 
jPapovich et al. 2004| [Rodighic ro et al. 2005| and at 850/im 
with SCUBA [e.g. Chapman et al. Il200;?r iresulting in the 



construction of much more accurate far-IR luminosity func- 
tions out to these redshifts with large number statistics. 

5.2 Expected Cosmic Far-Infrared Background 

The measurements of the infrared background have been 
carried out in several wavelength bands. Also, there has been 
rapid progress in resolving a significant fraction of the back- 
ground with deep galaxy counts at infrared wavelengths. 

The flux from extragalactic sources below the detec- 
tion limits create fluctuations in the background. The Cos- 
mic Far-Infrared Background (CFIRB) intensity Jcfirb pro- 
duced by all sources with the flux below the limiting flux 
Slim, is obtained from: 

In addition, the CFIRB fluctuations Pcfirb from sources 
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(b) Detection limits for Spitzer mission 

Figure 7. Comparison of detection limits with integration time for (a) AKARI mission and (b) Spitzer mission. The left panels are for 
SW band, the right panels for LW band. In the case of AKARI mission, the solid line shows the detection limits over an integration 
time for normal sampling mode and the dashed line shows those for sampling mode in high background region to avoid saturation. We 
plot two evolutionary models. Upper thick line is for the burst evolution model and lower thin line for the luminosity evolution model. 
Note that the integration time for single scan is 0.14s and 0.23s for the SW band and the LW band, respectively for AKARI and the 
integration time is plotted in logarithmic and linear scales for the AKARI and Spitzer missions, respectively. 



with a random distribution below a given detection limit 
Slim can be estimated as follows, 

PCFIRB = ^ ^ dS. (15) 

The detection limits Siim can be found in Table El 

We list the expected CFIRB intensity and resolved frac- 
tion for each mission in TableEl Dole et al. (|2004b|l predicted 
that the Spitzer mission can resolve 18% of the CFIRB at 
160 /im from their source count model. According to our es- 
timation with our source count model, the resolved CFIRB 
is expected to be 9 - 17% of the total CFIRB. As shown in 
Table|7| the Herschel & SPICA missions will resolve a much 
larger fraction of the CFIRB, i.e. 90 - 94% in the SW band 
and 60 - 72% in the LW band, compared to other missions. 

Table |H| lists the expected CFIRB fluctuations for each 
mission. We can not find a significant difference in the fluc- 
tuation levels of the Spitzer and the AKARI missions in the 
LW band since the detection limits already approach a flux 
range with a monotonic increase of fluctuation (see Fig. |3] 
for the burst evolution model). In the far-IR range, Lagache 
& Puget (I2000bll and Matsuhara et al. (2000) have stud- 
ied the detection of the CFIRB fluctuation with ISO data 
in Marano 1 region and the Lockman Hole, respectively. For 
comparison, we list the fluctuation estimated from our model 
for the ISO mission in Table |U] We find that the observed 
fluctuations are mostly located between the results from our 
two evolution models except for the observations at 90 fim. 
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Figure 9. Expected CFIRB fluctuations for the burst evolution 
model. Thick lines are for the LW band and thin lines for the SW 
band. Two arrows show the expected CFIRB fluctuation corre- 
sponding to the detection limits. 



Comparing with the model of Lagache, Dole & Puget 1200311 . 
our estimated fluctuations from the burst evolution model 
are in good agreement with those from their model. 

The power spectrum of cirrus emission at high Galactic 
latitudes (> 80 degree) also has a fluctuation of about 10® 
Jy^/sr at 0.01 arcmin"^ at 160 fim with a power index of 
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Figure 8. Expected optimal redshift distributions for the (a) Spitzer, (b) AKARI and (c) Herschel & SPICA missions. The 3 evolutionary 
models are shown for the SW and LW bands respectively at the confusion limit. This are the optimal distributions for each mission at 
the absolute sensitivity defined by our analysis. 



-2.9 ± 0.5. In order to distinguish the CFIRB fluctuations 
from the observed data that include both CFIRB and cirrus 
emission effectively, we require an area larger than 115, 125, 
and 19 arcmin^ for the Spitzer, AKARI, and the Herschel 
& SPICA missions, respectively for the case of luminosity 
evolution model. If we adopt the burst evolution model, the 
minimum required area for the resolution of the estimated 
power spectrum should increase to 463, 490, 36 arcmin^ for 
the Spitzer, AKARI, and the Herschel & SPICA missions. 



respectively. The shaded area in Figure ITUl covers the fluc- 
tuations for both evolutionary models. In this estimation, 
we do not consider the clustering of extragalactic sources. 
As sky confusion noise increases in the low Galactic lati- 
tude regions, detection limits should also degrade, and the 
CFIRB fluctuation becomes larger (see Figure \m . As can 
be seen in Figure [TTI the fluctuations show a monotonic in- 
crease to medium cirrus regions. Since the fluctuation power 
in low Galactic latitude regions should be more than 10^° 
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Table 7. Expected intensity, fluctuation, and resolved fraction of CFIRB for each mission. The upper flux is set to be the final confusion 
limit. 
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Intensity 
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Resolved 


Intensity 


Resolved 


Intensity 


Resolved 


Space Mission 


(MJy/sr) 


(%) 


(MJy/sr) 


(%) 


(MJy/sr) 


(%) 


(MJy/sr) 


(%) 


Spitzer 


0.063 


48 


0.33 


17 


0.12 


52 


0.69 


9 


AKARI 


0.094 


32 


0.31 


14 


0.21 


27 


0.66 


7 


Herschel & SPICA 


0.013 


90 


0.16 


60 


0.015 


94 


0.21 


72 



Table 9. Comparison of CFIRB fluctuations for ISO mission. 



A 


e 


'S'max 


Observations 


Predicted " 


Predicted (this work) 


(Atm) 


(arcmin) 


(mjy) 


(Jy2/sr) 


(JyVsr) 


(JyVsr) 


90 


0.4 - 20 


150 


13000 ± 3000 " 


5300 


2100 - 7200 


170 


0.6 - 4 


100 


7400 


12000 


3800 - 14000 


170 


0.6 - 20 


250 


12000 ± 2000 " 


16000 


5500 - 18000 



" Model from Lagache, Dole & Puget HD&H . 

^ Lower limit is estimated from luminosity evolution model and upper limit from burst evolution model. 
Observational analysis from Matsuhara et al. (20001. 
Observational analysis from Lagache &; Puget L2000bl . 



Table 8. Expected CFIRB fluctuations for each mission. 





Luminosity evolution 


Burst evolution 




(JyVsr) 


(JyVsr) 


Space Mission 


SW LW 


SW LW 



Spitzer 


35 


1600 


290 


12000 


AKARI 


130 


1800 


1200 


13000 


Herschel & SPICA 


0.4 


120 


0.6 


290 



Jy^/sr at 0.01 arcmin at 160 fim, it is difficult to extract 
the CFIRB fluctuations from the analysis of power spec- 
trum [Note that this CFIRB fluctuation is expected to have 
spatial structure below 1 arcmin"^ (see Figure ^J]. How- 
ever, we expect to detect CFIRB fluctuations in most of the 
low-to-medium cirrus regions. 



6 SUMMARY AND CONCLUSIONS 

In order to probe the confusion limit for infrared observa- 
tions, we generated source catalogues assuming a concor- 
dance (i.e. flat, dark energy dominated) cosmological world 
model {Ho = 72, Q = 0.3, A — 0.7) for 2 evolutionary sce- 



narios defined as the luminosity evolution model and burst 
model following Pearson II1996II and Pearson et al. II2001II . 
respectively. We also considered the sky confusion due to the 
infrared cirrus. Though the sky confusion is not a dominant 
noise source at high Galactic latitudes, we should take into 
account the effects of fluctuations of the sky brightness for 
large area surveys. 

Based upon the fluctuation analysis and photometry 
on simulated images, we find a composite estimator that 
represents the source confusion well. From our analysis of 
source confusion and source distribution models including 
galaxy evolution, we have estimated final confusion limits 
of 2.0 - 6.6 mJy and 25 - 67 mJy at 70 fim and 160 /im 
for the Spitzer mission, 6.4 - 17 mJy and 34 - 88 mJy 
at 75 ^m and 140 ^m for the AKARI mission, 0.077 - 
0.10 mJy and 2.3 - 4.4 mJy at 70 fim and 160 fim for 
the SPICA mission, and 0.077 - 0.10 mJy and 2.6 - 5.1 
mJy at 70 fim and 175 fim for the Herschel mission in low 
cirrus regions. If the source distribution follows the evolu- 
tionary models, the current and planned infrared missions 
will be mostly limited by source confusion. Other compo- 
nents affecting confusion are the fiuctuation from the zo- 
diacal hght, asteroids, and clustering of sources. As yet 
the fluctuation of the zodiacal hght on small scales has 
not been observed although we may expect this informa- 
tion from the current or next generation of space missions. 
Based upon the number distribution lITedesco et al. 2005t 
and SED ( |Miiller fc Lagerros 19981 IKim et al. 2003t of as- 
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Figure 10. Comparison between expected CFIRB and cirrus 
fluctuation for eacfi mission in tlie LW band. Tlie solid line shows 
the power spectrum of cirrus emission at high Galactic latitude. 
For comparison, we also plot the power spectrum in medium and 
high cirrus regions. The shaded area shows fluctuations and cor- 
responding spatial frequencies covering two evolution models at 
high Galactic latitude for each mission. The lower limit is for 
luminosity evolution and the upper limit for burst evolution. 



we have generated all sky maps for final confusion limits. 
These maps are shown in Appendix IHl 

We iiave also obtained the optimal (confusion limited) 
redshift distribution from each source count model. The 
redshift distributions for each mission show significant dif- 
ferences between the no evolution scenario and evolution- 
ary models. However, in order to distinguish between the 
two evolution scenarios, higher sensitivity /resolution mis- 
sions (e.g., Herschel and SPICA) are required compared to 
the present, relatively low resolution missions. We have es- 
timated that the CFIRB will be resolved to 48 - 52% and 
9 - 17% at 70 /Ltm and 160 ^m for the Spitzer mission, 27 
- 32% and 7 - 14% at 75 ^m and 140 /im for the AKARI 
mission, and 90 - 94% and 60 - 72% at 70 ^im and 160 ^.m 
for the Herschel & SPICA missions. We also found that we 
can detect the CFIRB fluctuations in most low-to-medium 
cirrus regions. 
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Figure 11. Expected CFIRB fluctuation as a function of cirrus 
mean brightness in the LW band. The two lines plotted for each 
mission are for the burst evolution model (upper line) and for the 
luminosity evolution model (lower line). 

teroids, we roughly estimated the effect of source confusion 
from asteroids. Since the slope of the flux-number distribu- 
tion is shallower than that of Euclidean space and the num- 
ber density at any given flux is much smaller than that of our 
evolutionary source distributions, we expect that the effect 
from asteroids near the ecliptic plane is not severe, i.e., less 
than 2% and 0.5% in the SW and LW bands, respectively. 
The fluctuation from the clustering of sources results in a 
10 ~ 15% increase in the confusion fTak eucIn fc Ishii 20041 
[Negrello et al. 2004| l , and can be comparable to the sky con- 
fusion in regions of medium cirrus brightness in the LW 
band. This effect of clustering on the confusion should be 
detectable in the SW band or in regions of low cirrus bright- 
ness. Using our study for sky confusion and source confusion. 
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APPENDIX A: PHOTOMETRIC RESULTS FOR 
SIMULATED DATA 

We have carried out photometry on our simulated images 
under various conditions. The results of the photometry is 
shown below. 

Fig. lAll shows the photometric results for each mission 
considering only source confusion. The results for the no evo- 
lution model are clearly segregated from those of the evolu- 
tion models due to the low source confusion. We also found 
that the completeness and the reliability fall rapidly since 
the source confusion is more severe for relatively low reso- 
lution observations such as is the case for the LW bands of 
AKARI and Spitzer. 

Fig. IA2I shows the photometric results for the AKARI 
mission with the no evolution model and various levels of 
cirrus background. Since the contribution from the sky con- 
fusion in the LW band is comparable to the source confusion 
(see the top- left panel of Fig.jSJ, the completeness and relia- 
bility are well segregated as the cirrus background increases. 

Fig. IA3I shows the photometric results for the AKARI 
mission for the two evolutionary models and various levels 
of cirrus background. Since the contribution from the source 
confusion in the LW band is stronger than that in no evolu- 
tion model (Fig. IA2l we find that the plots of completeness 
and reliability gather closer compared to those in the no 
evolution model. 

We also show the photometric results for the Spitzer 
mission in the case of the two evolutionary models in Fig. 
IA4I Although the cirrus background becomes higher, the 
completeness and reliability show similar values in the case 
of both the SW and LW bands with a mean brightness (B\) 
< 10 MJy sr~^ due to the relatively more severe source 
confusion compared to the case of the no evolution model. 
In addition, we found that though the AKARI and Spitzer 
missions have similar apertures (0.69m and 0.85m, respec- 
tively), the larger pixel size of AKARI makes both the sky 
and source confusion more severe. 

Fig. IA5I shows the photometric results for the Her- 
schel and SPICA missions from the two evolutionary mod- 
els. Since the source confusion is much more significant in 
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Figure Al. Photometric results for eacli mission considering only source confusion. Upper panels show the completeness and lower 
panels the reliability. Left and right panels show the SW band and LW band, respectively. A lighter and thinner line for the same line 
style means the result for a different source distribution evolutionary model. 
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Figure A3. Photometric results for the AKARI mission with the two evolutionary models. Upper panels show the completeness for 
various cirrus backgrounds and lower panels, the reliability. Left and right two panels show the results for the luminosity evolution and 
burst evolution models, respectively. Each line represents the photometric result for each cirrus background. Due to the equivalent or 
excessive sky confusion compared to the source confusion, each result shows a clear difference. 




Figure A4. Photometric results for the Spitzer mission from the two evolutionary models. The explanations for each panel and lines 
are the same as in Fig. IA3i Compared to the case of the luminosity evolution model, the more severe source confusion makes each result 
similar in the case of burst evolution model especially in the LW band. 




APPENDIX B: ALL SKY CONFUSION MAPS 

Using our results for sky confusion and source confusion, we 
have generated all sky maps for the final confusion limits. 
Since the final confusion limits in the SW band are domi- 
nated by source confusion, we show our results only for the 
case of the LW band (see Figs lBlllB2l and lB4ll except for all 
sky survey of the AKARI mission. Since the instrumental 
noise of the all sky survey of the AKARI mission decreases 
in proportion to the square root of the number of scans, we 
show here the final detection limit map including instrumen- 
tal noise by considering the visibility for the all sky survey'^ 
(see Figure IHst . 
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Figure A2. Photometric results for the AKARI mission with 
the no evolution model. Upper panels show the completeness for 
various cirrus backgrounds and lower panels, the reliability. Left 
and right panels show the results for the SW band and LW band, 
respectively. Each line represents the photometric result for each 
cirrus background. Due to the equivalent or excessive sky confu- 
sion compared to the source confusion, each result shows a clear 
difference. 



most of the cirrus background regions (see the lower panels 
of Fig. ISJ, the completeness does not change significantly, 
irrespective of the amount of cirrus background. However, 
from the reliability plot we see that the sky confusion still 
affects the measurement of the source flux. 



^ Further information can be found at the following url: 
http :/ / www. ir. isas.jaxa.jp / A KA RI/ Observation/vis / 




Figure A5. Photometric results for the Herschel & SPICA missions from the two evolutionary models. The explanations for each panel 
and lines are the same as in Fig. EH Since the source confusion is the dominant noise source, each result shows no significant difference, 
regardless of the cirrus background except for the case of the reliability in the LW band. 
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Figure Bl. Final confusion limit map for the Spitzer mission in Galactic coordinates with the Galactic center located at the center of 
image. The left and right panels show the results for the luminosity evolution model and the burst evolution model, respectively. 
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Figure B2. Final confusion limit map for the scan mode observation of the AKARI mission in Galactic coordinates with the Galactic 
center located at the center of image. The left and right panels show the results for the luminosity evolution model and the burst evolution 
model, respectively. 
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Figure B3. Final detection limit map for the survey mode observation of the AKARI mission in Galactic coordinates with the Galactic 
center located at the center of image. The left and right panels show the results for the SW band and the LW band, respectively. Since 
the detection limits mainly depend on the instrumental noise except for near the ecliptic poles and Galactic plane, we show only the 
case for the luminosity evolution model. 





Figure B4. Final confusion limit map for the Herschel & SPICA missions in Galactic coordinates with the Galactic center located 
at the center of image. The left and right panels show the results for the luminosity evolution model and the burst evolution model, 
respectively. Most of the regions are limited by source confusion. 



